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SUMMARY 

Several counterdiffusion type crystal growth experiments have been con- 
ducted In space. Improvements In crystal size and quality have been attributed 
to reduced natural convection In the microgravity environment. One series of 
experiments called DMOS (Diffusive Mixing of Organic Solutions) was designed 
and conducted by researchers at the 3M Corporation and flown by NASA on the 
space shuttle. Since only limited Information about the mixing process Is 
available from the space experiments, a series of ground based experiments 
were conducted to further Investigate the fluid flow within the DMOS crystal 
growth cell. Solutions with density differences In the range of 1 0—7 to 
10-4 g/cm3 were used to simulate microgravity conditions. The small density 
differences were obtained by mixing D 2 O and H 2 0. Methylene blue dye was 
used to provide flow visualization. The extent of mixing was measured photo- 
metrically using the 662 nm absorbance peak of the dye. Results Indicate that 
extensive mixing by natural convection can occur even under microgravity con- 
ditions. This Is qualitatively consistent with results of a simple scaling 
analysis. Quantitative results are In close agreement with ongoing computa- 
tional modeling analysis. 


1. INTRODUCTION 

Organic crystals have shown considerable potential for use In nonlinear 
optical devices (ref. 1). The usefulness of these materials In devices Is 
often limited by the Inability to grow large crystals with a sufficiently low 
density of defects. On Earth, temperature and composition gradients result In 
natural (bouyancy driven) convection, which can Interfere with crystal growth 
(ref. 2). One technique for growing crystals Involves the counterdiffusion of 
two reactants (A and B) and precipitation of the salt (AB) In the Interfaclal 
region. This approach Is particularly useful when the salt Is not sufficiently 
soluble In common solvents to allow growth by precipitation from a homogeneous 
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solution. In the counterdiffusion process natural convection can occur along 
with molecular diffusion as a result of density gradients within the container. 
The density gradients are a result of: (1) concentration gradients for com- 

ponents A and B, (2) the density difference between the solution and the solid 
salt crystal, and (3) thermal gradients resulting from heat released during 
salt formation and crystallization. These density gradients cannot be elimi- 
nated In practical crystal growth experiments. However, the natural convection 
caused by the density gradients can be considerably reduced If the crystal 
growth experiment Is conducted In space. 

Organic crystals were grown by researchers at the 3M Corporation In a 
series of space experiments called DMOS (Diffusive Mixing of Organic Solu- 
tions). Included In DMOS-2 were two experiments designed to measure the 
extent of solution mixing In microgravity. In order to examine the fluid 
dynamics more closely, a series of ground based experiments were designed to 
simulate the fluid flow which occurred during the DMOS experiments. These 
experiments were conducted at the NASA Lewis Research Center In collaboration 
with researchers at the 3M Corporation and the NASA Marshall Space Flight 
Center. This paper describes the details of the ground based experiments and 
results which apply both to the DMOS experiments and to other similar experi- 
ments with horizontal composition gradients. 


2. EXPERIMENTAL PROCEDURES 

The DMOS crystal growth apparatus consists of three cylindrical chambers 
placed end to end. Between the chambers a sliding gate valve assembly Is used 
to control mixing. For the ground based experiments described In this paper, 
the simpler two chamber cell shown In figure 1 was used. The chambers of this 
cell are rectangular for better photographic flow visualization, but the gate 
valve assembly Is Identical to that used In the DMOS experiments. In a typical 
experiment the left chamber Is filled with a 0.01 g/1 solution of methylene 
blue dye In water while the right chamber Is filled with pure water. The den- 
sity difference between fluids In the chambers Is controlled by adding deuter- 
ated water (D 2 O) to one chamber or the other. Qualitative Information about 
the flow patterns during mixing was obtained using time lapse photography at 2 
frames per second. Quantitative results were obtained by photometric analysis 
using the 662 nm absorbance peak of the methylene blue dye. This was done by 
closing the gate valve after the desired mixing time, shaking the apparatus to 
get a uniform dye concentration, and taking samples from each chamber for 
analysis. At sufficiently long times the mass fraction of dye In each chamber 
Is 0.5. If f(t) Is the mass fraction of dye In the right chamber at any 
time, t, the extent of mixing, Xm(t), Is defined by equation 1. 


%Xm = 100X 



( 1 ) 


Using this definition, the extent of mixing goes from 0 to 100 percent as the 
mixing time goes from zero to Infinity. The temperature was maintained at 
30±0.1 °C during the the mixing experiments by Immersing the apparatus In a 
circulating water bath. A quartz thermometer was used to measure the tempera- 
ture difference between various points In the bath. When the density differ- 
ence between the solutions was less than 10-5 g/cm3, heating of the dye solu- 
tion by the background lighting used for photography began to Influence the 
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fluid flow. Therefore, when quantitative data was taken, no background light- 
ing was usfed, and the water bath was covered with a black cloth to exclude 
room light. 

The diffusion coefficient of the methylene blue dye In O 2 O was measured 
by Prof. Ernst von Heerwall at the University of Akron using a pulsed gradient 
spin echo (PGSE), nuclear magnetic resonance (NHR) technique. Physical prop- 
erties of H 2 O and D 2 O used In section 3 were obtained from standard tables 
using values at 30 ®C. 


3. RESULTS AND DISCUSSION 
3.1 The Microgravity Environment 

Experiments conducted on the shuttle experience both translational and 
rotational accelerations from many sources (ref. 3). Atmospheric drag causes 
a steady acceleration of the center of mass antiparallel to the orbit velocity 
A steady acceleration perpendicular to the orbit also exists for points not 
located on the orbital path of the center of mass. Oscillatory accelerations 
occur as a result of structural vibrations, and transient accelerations occur 
as a result of attitude maneuvers, satellite deployments, crew motions, etc. 

An experiment can be run during quiet periods, so that only steady accelera- 
tions are Important. In this paper only flows resulting from a steady accel- 
eration perpendicular to the Initial density gradient are considered. 


3.2 Microgravity Simulation-Experiments 

It will be assumed for now that the fluid flow which occurs during a 
space experiment can be simulated In a ground-based experiment If the product 
of the solution density difference and the effective gravitational accelera- 
tion Is the same In each experiment. The validity of this assumption will be 
oiscusse^ wrm p r^ detan In section 3.3. A detailed description of the accel- 
erations which occurred during the DMOS-2 experiments will be published later. 
Preliminary analysis of accelerometer data Indicates that the effective accel- 
eration was on the order of 10-5 times that on Earth. Density differences 
In ground-based simulations must therefore be smaller than those In the space 
experiments by the same factor. The solutions used to obtain these small 
density differences must be different from the solutions used In the space 
experiments. In some cases solutions very similar to those used In the space 
experiment can be used In ground-based simulations If Isotopic mixtures (such 
as D 2 O/H 2 O) are used to control solution densities. This allows the density 

difference between two solutions to be varied over many orders of magnitude 
with only a minor modification of chemical properties. 

An experiment was conducted In the apparatus of figure 1 using a blue 
solution (H 2 O + dye) In the left chamber and pure water (H 2 O) In the right 

chamber. A dye concentration of 0.01 g/1 was needed for photographic flow 
visualization. When the gate valve was opened, the blue solution flowed from 
left to right through the three bottom channels (counterclockwise flow). This 
Indicates that the dye causes an Increase In the density of the water. How- 
ever, the density change Is too small to be measured directly by standard 
techniques. In order to determine the effect of the dye on the density of the 
water, a series of experiments were run In which the density of the solution 
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In the right chamber was Incrementally Increased by adding various amounts of 
O 2 O to thfe water. Each experiment was stopped after 2 hr, and the extent 
of mixing was determined. Results of these experiments are shown In figure 2. 
Positive and negative values Indicate flow In the counterclockwise and clock- 
wise directions respectively. If Ideal mixing Is assumed, the density differ- 
ence between the H 2 O/D 2 O solution and pure H 2 O can be calculated using 
equation 2. 


Ap = V D 2 0^D 2 0 ' p H 2 oj 

v„ „ volume fraction of DO 
0 2 0 2 

p h Q density of pure H 2 0 

Q density of pure D 2 0 


( 2 ) 


According to figure 2, the effect of the dye Is equivalent to a D 2 O 
volume fraction of 3.5x10-5. Using this value for v D2 g and (pq 2 q " ^ 0 ) = 
0.108 g/cm3 in equation 2, the density Increase caused by adding 0.01 g/1 of 
dye to H 2 0 Is calculated to be 3. 8±0. 1x10-6 g/cm3. Notice that It Is possi- 
ble to prepare solutions with density differences on the order of 10-6 g/cm3 
without knowing the absolute densities of the solutions to the same accuracy. 
This Is possible only If the same source of water Is used to prepare both 
solutions and no significant concentration of Impurities Is Introduced during 
preparation of the solutions. The consistency and reproducibility of results 
obtained throughout this work suggest that these conditions were satisfied. 
Another potential source of error Is thermally Induced density gradients. The 
maximum temperature difference measured between any two points In the water 
bath Just outside of the apparatus was 0.007 °C. Using 3.0x10-4 K-l as the 
volumetric thermal expansion coefficient of water, the maximum possible ther- 
mally Induced density difference Is 2.1x10-6 g/cm3. Since the time averaged 
temperature difference between any two points outside of the apparatus Is much 
less than 0.007 °C, thermally Induced density differences Inside the apparatus 
are probably on the order of 10-7 g/cm3 or less. 


In order to observe qualitatively how mixing occurs, a series of experi- 
ments with density differences ranging from 10-6 to 10-4 g/cm3 were conducted. 
The flow was qualitatively the same throughout this range of density differ- 
ences. The denser (blue) solution flowed from left to right through the bot- 
tom three channels of the gate valve assembly. When the blue solution entered 
the right chamber. It flowed like a waterfall In a layer about 0.2 cm thick 
along the vertical wall of the gate valve assembly to the bottom of the cham- 
ber. The flow continued until the levels of the blue solution were nearly 
equal In each chamber. 

In some of the DMOS space experiments the density difference between 
solutions In different chambers was on the order of 1 g/cm3. in order to 
simulate the flow which occurred In these space experiments, a series of 
ground based experiments with Ap = 4.9x10-6 g/cm3 were run. The extent of 
mixing as a function of time Is shown In figure 3. Nearly complete mixing Is 
obtained In about 32 hr. At this time each chamber contains the blue solution 
In the bottom of the chamber and a clear solution at the top of the chamber 
with a diffuse horizontal Interface separating the two solutions. No observa- 
ble change In the dye distribution was observed between 32 and 56 hr. 
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3.3 Microgravity Simulation-Theory 


In order to describe the fluid flow which occurs during experiments con- 
ducted In the apparatus of figure 1, It Is convenient to consider flow through 
the valves and flow within the chambers separately. As a simplification, flow 
through the valves can be described as Polseullle flow through a rectangular 
channel. The average velocity of fluid within the valves can then be calcu- 
lated from equation 3 (ref. 4). 


AP 

L 




(3) 


u average fluid velocity 

AP pressure difference 
L channel length 


p fluid density 

d hydraulic diameter 

X friction factor 


The hydraulic diameter Is defined as 4 x (area)/(c1rcumference) for the 
rectangular channels. Using the Reynolds (Re = ud/v) number based on the 
hydraulic diameter, the friction factor Is approximately X = 64/Re. If the 
hydrostatic pressure difference resulting from a density difference of 
4.9x10-6 g/cm3 Is used In equation 3, the fluid velocities through the bottom 
three channels are 3.02x10-3, 2.14x10-3, and 0.73x10-3 cm/s respectively, 
beginning with the bottom channel. For these velocities the entrance length 
for fully developed flow [(He = (0.035)(d)(Re)] Is much less than the channel 
length. This condition and an assumption of quasisteady flow are required In 
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left to the right chamber Is calculated using these flow velocities, the extent 
of mixing after 3162 sec (log t = 3.5) Is calculated to be 18.7 percent. This 
Is close to the experimental value of 12.0 percent In figure 3. Numerical 
computations and computer models are currently being developed to provide a 
more accurate description of the flow. 


Since the denser solution emerging from the channels flows like a water- 
fall to the bottom of the chamber, the vertical flow velocity must be much 
larger than the horizontal velocity. This can be explained using simple 
scaling arguments (ref. 5). The relative Importance of various transport 
processes Is determined by calculating the magnitude of appropriate dimension- 
less parameters. The parameters of Interest are the Grashof number, Gr, and 
the Schmidt number. Sc, which are defined below. 

, x. 3 

Gr = 9(V)L (4) 

pv 

Sc = * (5) 
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g gravl Rational acceleration 
p average density 

Ap density difference 

L characteristic length 

v kinematic viscosity 

0 diffusion coefficient of dye 

Using Ap * 4.9x10-6 g/cm3, p * 1.00 g/cm3, g =» 980 cm/s2, L * 4.45 cm, 
v * 0.00804 cm2/s, and D * 3x10-6 cm2/s In equations 4 and 5 gives, Gr * 
6.6x103, Sc = 2.7x103, and ScGr = 1.8x107. An order of magnitude estimate of 
the vertical fluid velocity can be calculated from equation 6 or 7 (ref. 5). 

v = Gr1/2 (f) (for Gr >> (6) 

v = Gr^ (for Gr « 1) (7) 

Using equation 6 the vertical fluid velocity Is 0.146 cm/s. As expected, 
this Is much larger than the horizontal velocities calculated above. A bound- 
ary layer thickness (6) of 0.75 cm Is estimated using equation 8 (ref. 5). 


i ] 

L = Gr 1/4 


( 8 ) 


This Is within order of magnitude agreement with the observed thickness of 
0.2 cm. 


So far It has been Implicitly assumed that the motion of the dye molecules 
represents the motion of the fluid. The validity of this assumption depends on 
ScGr as described In equations 9 and 10. 


ScGr « 1 

/diffusion! 

(9) 


ydomlnant J 


ScGr » 1 

/ convect1on\ 

(10) 


1 dominant / 



Since ScGr » 1, transport of the dye by molecular diffusion Is much 
slower than by convection, and the motion of the dye does represent convection 
of the fluid. Diffusion Is not Important unless (Ap)g « 2.7x10-10 0/cm2-s2. 
This requires matching densities to within 10 -7 g/cm3 for a space experiment or 
10-13 g/cm3 for a ground based experiment. Use of such small density differ- 
ences Is not experimentally feasible. Furthermore, the time scale of a diffu- 
sion controlled experiment would be too long to be of practical Interest. 
Therefore, changes In the geometry of the apparatus In figure 1 would be 
required In order to conduct a diffusion controlled experiment. 

It remains to be determined whether or not matching the (Ap)g product 
Is sufficient for simulation of fluid flows, within the chambers It has been 
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shown that molecular diffusion Is not Important, so that the only Important 
dlmenslonlfess group Is Gr. For Polseullle flow between the chambers the mass 
transfer rate Is proportional to (Ap)g/pw. Since flow between chambers 
depends on viscosity as l/t> and flow within the chambers depends on viscosity 
as l/t>2, the overall process can not be simulated unless the solution viscosi- 
ties are the same In both ground based and space experiments. However, both 
the viscosity and density of most organic and aqueous solutions containing 
moderate concentrations of low molecular weight solutes do not vary by more 
than one order of magnitude. As a result, matching the (Ap)g product should 
provide reasonable (but not exact) simulation In most cases of Interest. 


CONCLUSIONS 

The fluid flow which occurs during counterdiffusion type crystal growth 
experiments In space can be studied using ground based simulations. For Iso- 
thermal mixing of solutions, approximate simulation can be accomplished by 
matching the (Ap)g product. Density differences as small as 10-6 g/cm3 can 
be easily obtained using Isotopic mixtures. For experiments with a geometry 
similar to that In figure 1, natural convection Is expected to be the dominant 
mode of transport for both space and ground based experiments If the direction 
of the gravitational acceleration Is perpendicular to that of the Initial den- 
sity gradient. The experimental approach to ground based simulation and the 
scaling analysis described In this paper can be generally applied to other 
systems with horizontal composition gradients In a vertical gravitational 
field. 
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FIGURE 1.- CROSS SECTION OF 2 CHAMBER MIXING CELL 
USED FOR GROUND-BASED EXPERIMENTS. 
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FIGURE 2.- EXTENT OF MIXING AFTER TWO 
HOURS VERSUS VOLUME FRACTION OF D 2 0 
IN RIGHT CHAMBER. LEFT CHAMBER CON- 
TAINS 0.01 G/l METHYLENE BLUE DYE. 
h 2 o IS THE SOLVENT IN BOTH CHAMBERS. 




FIGURE 3.- EXTENT OF MIXING VERSUS 
TIME FOR INITIAL DENSITY DIFFERENCE 
BETWEEN CHAMBERS OF 4.9x10~ 6 g/CM 3 . 
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